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Abstract 

Microwave and submillimeter molecular transition frequencies between nearly degenerated rota- 
tional levels, tunneling transitions, and mixed tunneling-rotational transitions show an extremely 
high sensitivity to the values of the fine-structure constant, a, and the electron-to-proton mass ra- 
tio, [i. This review summarizes the theoretical background on quantum-mechanical calculations of 
the sensitivity coefficients of such transitions to tiny changes in a and \i for a number of molecules 
which are usually observed in Galactic and extragalactic sources, and discusses the possibility 
of testing the space- and time-invariance of fundamental constants through comparison between 

^O precise laboratory measurements of the molecular rest frequencies and their astronomical counter- 

parts. In particular, diatomic radicals CH, OH, NH + , and a linear polyatomic radical C3H in II 
electronic ground state, polyatomic molecules NH3, ND3, NH2D, NHD2, H2O2, H30 + , CH3OH, 
and CH3NH2 in their tunneling and tunneling-rotational modes are considered. It is shown that 
sensitivity coefficients strongly depend on the quantum numbers of the corresponding transitions. 
This can be used for astrophysical tests of Einstein's Equivalence Principle all over the Universe at 

-t— > an unprecedented level of sensitivity of ~ 10 -9 , which is a limit three to two orders of magnitude 

lower as compared to the current constraints on cosmological variations of a and fi: Aa/a < 10~ 6 , 
A/j,/n < 1CT 7 . 
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I. INTRODUCTION 

The fundamental laws of particle physics, in our current understanding, depend on 28 
constants including the gravitational constant, G, the mass, m e , and charge, e, of the elec- 
tron, the masses of six quarks, m u , m<j, m c , m 8 , m t , and m^, the Planck constant, h, the 
Sommerfeld constant a, the coupling constants of the weak, g w , and strong, g s , interactions, 
etc. The numerical values of these constants are not calculated within the Standard Model 
and remain, as Feynman wrote about the fine structure constant a in 1985, "one of the 
greatest mysteries of physics" [lj. However, it is natural to ask whether these constants are 
really constants, or whether they vary with the age of the universe, or over astronomical 
distances. 

The idea that the fundamental constants may vary on the cosmological time scale has 
been discussing in different forms since 1937, when Milne and Dirac argued about possible 
variations of the Newton constant G during the lifetime of the universe [21 13|. Over the past 
few decades, there have been extensive searches for persuasive evidences of the variation 
of physical constants. So far, there was found no one of them. The current limits for 
dimensionless constants such as the fine structure constant, a = e 2 /hc, and the electron 
to proton mass ratio, jjl = m e /m p , obtained in laboratory experiments and from the Oklo 
natural reactor are on the order of one part in 10 15 — 10 17 J1H6] and one part in 10 14 — 10 16 [7]- 
[9] per year, respectively. The detailed discussion of ideas behind laboratory experiments can 
be found in a review |1UJ . 

Assuming that the constants are linearly dependent on the cosmic time, the same order 
of magnitude constraints on the fractional changes in Aa/a = (a bs — &i&b)/ a iah and in 
A/i//i = (/x bs — A*iab)//^iab are stemming from astronomical observations of extragalactic 
objects at redshifts z ~ 1 — 5 [l~lTfi"5] . Less stringent constraints at a percent level have been 
obtained from the cosmic microwave background (CMB) at z ~ 10 3 [T6HT8] and big bang 
nucleosynthesis (BBN) at z ~ 10 10 [191 120] • We note that space and/or time dependence 
of a based on optical spectra of quasars and discussed in the literature [211 an d references 
therein] is still controversial and probably caused by systematic effects since independent 
radio-astronomical observations, which are more sensitive, show only null results for both 
Aa/a and A///// [221123]. 

Surprisingly, it looks as if the Einstein heuristic principle of local position invariance 
(LPI) - - the outcome of any local non- gravitational experiment is independent of where and 
when in the universe it is performed --is valid all over the universe, i.e., at the level of 
~ 10~ 6 neither a no /i deviate from their terrestrial values for the passed 10 10 yr. In the 
Milky Way, it was also found no statistically significant deviations of Afi/fi from zero at 
even more deeper level of ~ 10~ 8 [24-26J. 

However, the violation of the LPI was predicted in some theoretical models such as, for 
example, the theory of superstrings which considers time variations of a, g w , and the QCD 
scale Aq CD (i.e., /x since m p ex Aq CD ) and thereby opening a new window on physics beyond 
the Standard Model [271 an d references therein]. If the fundamental constants are found 
to be changing in space and time, then they are not absolute but dynamical quantities 
which follow some deeper physical laws that have to be understood. Already present upper 
limits on the variation of the fundamental constants put very strong constraints on the 
theories beyond the Standard Model [28, and references therein]. This motivates the need 
for more precise laboratory and astronomical tests of the LPI. Of course, there are also other 
attempts to look for the new physics. For example the electric dipole moments (EDMs) of 



the elementary particles are very sensitive to the different extensions of the Standard Model. 
Present limit on the EDM of the electron significantly constrains supersymmetrical models 
and other theories |2Sl ED] • 

In this review we will consider tests of LPI which are based on the analysis of microwave 
and submillimeter 1 astronomical spectra and which are essentially more sensitive to small 
variations in a and [i than the test based on optical spectral observations of quasars. 

II. DIFFERENTIAL MEASUREMENTS OF Aa/a AND Afi/fi FROM ATOMIC 
AND MOLECULAR SPECTRA OF COSMIC OBJECTS 

Speaking about stable matter, as, for example, atoms and molecules, we have only seven 
physical constants that describe their spectra [31] : 

G, Aqcd, a, m e , m u , m d , m s . 

The QCD scale parameter Aqcd and the masses of the light quarks u, d, and s contribute 
to the nucleon mass m p (with Aqcd ^> itt-u + "^d + Wb) and, thus, the electron-to-proton 
mass ratio fi is a physical constant characterizing the strength of electroweak interaction in 
terms of the strong interaction. 

In the nonrelativistic limit and for an infinitely heavy pointlike nucleus all atomic tran- 
sition frequencies are proportional to the Rydberg constant, R, and the ratios of atomic 
frequencies do not depend on any fundamental constants. Relativistic effects cause correc- 
tions to atomic energy, which can be expanded in powers of a 2 and a 2 Z 2 , the leading term 
being a 2 Z 2 R, where Z is atomic number. Corrections accounting for the finite nuclear mass 
are proportional to R^/Z, but for atoms they are much smaller than relativistic corrections. 

Astronomical differential measurements of the dimensionless constants a and /i are based 
on the comparison of the line centers in the absorption/emission spectra of cosmic objects 
and the corresponding laboratory values. It follows that the uncertainties of the laboratory 
rest frequencies and the line centers in astronomical spectra are the prime concern of such 
measurements. It is easy to estimate the natural bounds set by these uncertainties on the 
values of Aa/a and A/i//i. 

Consider the dependence of an atomic frequency u on a in the comoving reference frame 
of a distant object located at redshift z J32J 133] : 

u z = u + qx + 0(x 2 ), x = (a z /a) 2 — l. (1) 

Here u and u z are the frequencies corresponding to the present-day value of a and to a 

change a — > a z at a redshift z. In this relation, the so-called q factor is an individual 

parameter for each atomic transition. 

If a z 7^ a, the quantity x in ([I]) differs from zero and the corresponding frequency shift 

Acu = u z — uj is given by 

Aw Act 

= Q—, 2 

to a 

1 The frequency range 1 GHz < v < 300 GHz is usually referred to as a microwave range. Molecular 
transitions below 1 GHz (wavelength A > 30 cm) are from a low-frequency range which is restricted by 
the ionospheric cut-off at 10 MHz (A = 30 m). 



where Q = 2q/uj is the dimensionless sensitivity coefficient and Aa = (a z — a) /a is the 
fractional change in a. Here we assume that |Aa/a| <C 1. The condition a z ^ a leads to a 
change in the apparent redshift of the distant object Az = z — z: 



Aoj Az Av 



U) 1 + z 



(3) 



where Av is the Doppler radial velocity shift. 

If u/ is the observed frequency from the distant object, then the true redshift is given by 



!+< = £, (4) 



whereas the shifted (apparent) value is 

l + i=£. (5) 

Now, if we have two lines of the same element with the apparent redshifts Z\ and z% and the 
corresponding sensitivity coefficients Qi an d Qi, then 

Aa ^ -z 2 Av 

AQ = = — . 6 

a 1 + 2 c 

Here At> = V\ — t>2 is the difference of the measured radial velocities of these lines, and 
AQ = Q 2 — Qi is the corresponding difference between their sensitivity coefficients. By 
comparing the apparent redshifts of two lines with different sensitivity coefficients Q we can 
study variation of a on a cosmological timescale. 

Unfortunately, optical and UV transitions of atoms and molecules are not very sensitive 
to changes in a and /1. The sensitivity coefficients of atomic resonance transitions of usually 
observed in quasar spectra chemical elements (C, N, O, Na, Mg, Al, Si, S, Ca, Ti, Cr, 
Mn, Fe, Co, Ni, Zn) are very small, Q ~ {&Z) 2 <C 1 |34j . The same order of magnitude 
sensitivity coefficients to /i variations have been calculated for the UV transitions in the 
Lyman and Werner bands of molecular hydrogen H2 [35H37] , and for the UV transitions in 
the 4th positive band system A 1 !! — X 1 S + of carbon monoxide CO |38j. 

Small values of Q and AQ put tough constraints on optical methods to probe Aa/a and 
Afi/fi. Let us consider an example of Fell lines arising from the ground state 3<i 6 ( 5 .D)4s. 
In quasar spectra we observe 7 resonance transitions ranging from 1608 A to 2600 A with 
both signs sensitivity coefficients: Qai608 — —0.0322, Qai6h — +0.0502, and Q ~ +0.08 
for transitions with A > 2000 A [3H note a factor of two difference in the definition of the 
coefficients Q with the present work]. This gives us the maximum value of AQ ~ 0.11 which 
is known with an error of ~ 30%. From (pi) it follows that a variance of Aa/a ~ 10~ 5 
would induce a velocity offset At; ~ 0.3 km s _1 between the 1608 A line and any of the line 
with A > 2000 A. We may neglect uncertainties of the rest frame wavelengths since they 
are ~ 0.02 km s _1 [40J. If both iron line centers are measured in quasar spectra with the 
same error a v , then the error of the offset Av is a^ v = v^V The error a& v is a statistical 
estimate of the uncertainty of Av, and, hence, it should be less than the absolute value 
of Av. This gives us the following inequality to adjust parameters of spectral observations 
required to probe Aa/a at a given level: 

AQ Aa 
a v < c . (7) 

a/2 a 



At Aa/a ~ 10~ 5 , the required position accuracy should be a v < 0.25 km s _1 . A typical 
error of the line center of an unsaturated absorption line in quasar spectra is about l/10th 
of the pixel size (the wavelength interval between pixels) |41j . Current observations with the 
UV- Visual Echelle Spectrograph (UVES) at the ESO Very Large Telescope (VLT) provide 
a pixel size AA pix ~ 0.05 — 0.06 A, i.e., at A ~ 5000 A the expected error a v should be ~ 0.3 
km s _1 , which is comparable to the velocity offset due to a fractional change in a at the level 
of 10 -5 . Such a critical relationship between the 'signal' (expected velocity offset Av) and 
the error a v hampers measuring Aa/a at the level of ~ 10~ 5 from any absorption system 
taking into account all imperfections of the spectrograph and the data reduction procedure. 
Systematic errors exceeding 0.5 km s _1 are known to be typical for the wavelength calibration 
in both the VLT/UVES and Keck/HIRES spectrographs [121 HUH]. At this level of the 
systematic errors an estimate of Aa/a from any individual absorption-line system must be 
considered as an upper limit but not a 'signal'. Otherwise, a formal statistical analysis of 
such values may lead to unphysical results (examples can be found in the literature). 

The UV molecular spectra of H_2 and CO observed at high redshifts in the optical wave- 
length band encounter with similar difficulties and restrictions. The maximum difference 
between the sensitivity coefficients in case of H_2 is AQ ~ 0.06, the rest frame wavelength 
uncertainties are negligible, ~ 5 x 10~ 9 [15], and with the current spectral facilities at giant 
telescopes it is hard to get estimates of A/z//i at a level deeper than 10~ 5 . For carbon monox- 
ide such measurements have not been done so far but the expected limit on A/i//i should be 
> 10 -5 since CO lines are much weaker than H 2 [46J and therefore their line centers are less 
certain. The analogue of Eq.([6j) for the //-estimation from a pair of molecular lines is |47j : 

A /i _ Av v 1 - v 2 , . 

/i ~ cAQ ~ c(Q 2 - Qi) ' U 

and for a given level of A/x//x, molecular line centers should be measured with an error 

a v < -j=—c . (9) 

This means that at A///// ~ 10~ 5 , the required position accuracy should be a v < 0.13 km s _1 , 
or the pixel size AA P j x < 0.017 A at 4000 A. This requirement was realized in the VLT/UVES 
observations of the quasar Q0347-383 [48J where a limit on Afi/fi of (4.3 ± 7.2) x 10~ 6 was 
set. 

At present the only way to probe variation of the fundamental constants on the cosmo- 
logical timescale at a level deeper than 10~ 5 is to switch from optical to far infrared and 
microwave bands. In the microwave, or submillimeter range there are a good deal of molec- 
ular transitions arising in Galactic and extragalactic sources. Electronic, vibrational, and 
rotational energies in molecular spectra are scaled as E e \ : E v fo : E rot = 1 : /i 1 / 2 : jj. In other 
words, the sensitivity coefficients for pure vibrational and rotational transitions are equal to 
Q^ = 0.5 and Q^ = 1, respectively. Besides, molecules have fine and hyperfine structures, 
A-doubling, hindered rotation, accidental degeneracy between narrow close-lying levels of 
different types, which have a specific dependence on the physical constants. The advantage 
of radio observations is that some of these molecular transitions are approximately 100-1000 
times more sensitive to variations of ft and/or a than optical and UV transitions. 

In the far infrared waveband also lie atomic fine-structure transitions, which have sensi- 
tivity to a- variation Q a w 2 [19] . We can combine observations of these lines and rotational 



molecular transitions to probe a combination F = a 2 //i [50]. Besides, radio-astronomical 
observations allow us to measure emission lines from molecular clouds in the Milky Way with 
an extremely high spectral resolution (channel width ~ 0.02 km s _1 ) leading to stringent 
constraints at the level of ~ 10 -9 [23]. The level 10~ 9 is a natural limit for radio-astronomical 
observations since it requires the rest frequencies of molecular transitions to be known with 
an accuracy better than 100 Hz. At the moment only ammonia inversion transitions and 
18 cm OH A-doublet transitions have been measured in the laboratory with such a high 
accuracy [5T], [52] . 

In the next sections we consider in more detail the sensitivities of different types of 
molecular transitions to changes in a and ft. We are mainly dealing with molecular lines 
observed in microwave and submillimeter ranges in the interstellar medium, but a few low- 
frequency transitions with high sensitivities are also included in our analysis just to extend 
the list of possible targets for future studies at the next generation of large telescopes for 
low-frequency radio astronomy. 

III. DIATOMIC RADICALS IN THE II GROUND STATE: CH, OH, AND NH+ 

We start our analysis of the microwave spectra of molecules from the simplest systems 
- diatomic molecules with nonzero projection of the electronic angular momentum L on 
the molecular axis. Several such molecules are observed in the interstellar medium. Here we 
will mostly focus on the two most abundant species — CH and OH. Recently it was realized 
that A-doublet transitions in these molecules have high sensitivity to the variation of both a 
and fj, [53H55] . There are also several relatively low frequency transitions between rotational 
levels of the ground state doublet II1/2 and n 3 / 2 with sensitivities, which are significantly 
different from the typical rotational ones [56]. Then we will briefly discuss the NH + radical 2 , 
which is interesting because it has very low lying excited electronic state 4 E~. This leads 
to an additional enhancement of the dimensionless sensitivity coefficients Q [58J . The latter 
are defined as follows: 

Acu ^ Aa „ Aa , , 

— = Q a — + Q»— ■ 10 

uj a jj, 

A. A-doubling and ^-doubling 

Consider electronic state with nonzero projection A of the orbital angular momentum on 
the molecular axis. The spin-orbit interaction couples electron spin S to the molecular axis, 
its projection being E. To a first approximation the spin-orbit interaction is reduced to the 
form H so = AAH. Total electronic angular momentum J e = L + S has projection Q on the 
axis, Q = A + E. For a particular case of A = 1 and S = | we have two states H\/2 and n 3 / 2 
and the energy difference between them is: E(U 3 /2) — E(Ui/2) = A. 

Rotational energy of the molecule is described by the Hamiltonian: 

H rot = B(J-J e ) 2 (11a) 

= BJ 2 -2B(JJ e ) + BJ 2 , (lib) 

2 NH + has not yet been detected in space, its fractional abundance in star-forming regions is estimated 
7V(NH+)/7V(H 2 ) < 4 x 1CT 10 [57]. 



where B is the rotational constant and J is the total angular momentum of the molecule. 



The first term in expression (lib) describes conventional rotational spectrum. The last 
term is constant for a given electronic state and can be added to the electronic energy. 3 The 
second term describes fl-doubling and is known as the Coriolis interaction Hq ov . 



If we neglect the Coriolis interaction, the eigenvectors of Hamiltonian (11) have definite 
projections M and fl of the molecular angular momentum J on the laboratory axis and 
on the molecular axis respectively. In this approximation the states \J, M, A, E, fl) and 
| J, M, — A, — E, — fl) are degenerate, -Ej,±c = BJ(J + 1). The Coriolis interaction couples 
these states and removes degeneracy. New eigenstates are the states of definite parity p = 
±1 ESI: 



\j,M,n,p)= (\j,m,q) + P (-i) J - s \j,M,-n}) /V2. (12) 

The operator Hqot can only change quantum number fl by one, so the coupling of states 
|fl) and | — fl) takes place in the 2fl order of the perturbation theory in Hq ov . 

The fl-doubling for the state IIx/2 happens already in the first order in the Coriolis 
interaction, but has additional smallness from the spin-orbit mixing. The operator i^cor can 
not directly mix degenerate |A, E, fl) states |1, — |, |) and | — 1, |, — |) because it requires 
changing A by two. Therefore, we need to consider spin-orbit mixing of the II and E states: 

|n = i> = |i,-ii> + C|o,y>, (is) 

where 

C~A/(£ri-£ s ), ( 14 ) 

and then 

(fl = ii# Cor |fi = -|> = 2CB(J+ §)(A = l|L a |A = 0). (15) 

Note that ( depends on the non-diagonal matrix element (ME) of the spin-orbit interaction 



and Eq. (14) is only an order of magnitude estimate. It is important, though, that non- 
diagonal and diagonal MEs have similar dependence on fundamental constants. We conclude 
that fl-splitting for the Ili/ 2 level must scale as ABJ/(Eu — E^,). The fl-doubling for II3/2 
state takes place in the third order in the Coriolis interaction. Here i^cor has to mix first 



states II3/2 with II1/2 and H-3/2 with n_i/ 2 before ME (15) can be used. Therefore, the 
splitting scales as B 3 J 3 /[A(En — E^)]. 

The above consideration corresponds to the coupling case a, when \A\ ^> B. In the 
opposite limit the states Ili/2 and II3/2 are strongly mixed by the Coriolis interaction and 
spin S decouples from the molecular axis (coupling case b). As a result, the quantum 
numbers E and fl are not defined and we only have one quantum number A = ±1. The 
A-splitting takes place now in the second order in the Coriolis interaction via intermediate 
E states. The scaling here is obviously of the form B 2 J 2 /(En — E-%). Note that in contrast 
to the previous case \A\ ^> B, the splitting here is independent on A. 



Note that this term contributes to the separation between the states n^ and n 3 / 2 . This becomes 
particularly important for light molecules, where the constant A is small. 



We can now use found scalings of the A- and fi-doublings to determine sensitivity coef- 
ficients (10). We only need to recall that in atomic units A oc a 2 and Boc/i. We conclude 



that for the case a the ^-doubling spectrum has following sensitivity coefficients: 

State 2 II 1/2 : Q a = 2, Q ft = l, (16a) 

13/2 : Q a = — 2 , Q[i 



State 2 n 3/2 : Q a = -2 , Q M = 3 . (16b) 



For the case b, when S is completely decoupled from the axis, the A-doubling spectrum has 
following sensitivity coefficients: 

State n : Q a = , Q^ = 2 . (16c) 

When constant A is slightly larger than B, the spin S is coupled to the axis only for lower 
rotational levels. As rotational energy grows with J and becomes larger than the splitting 
between states II1/2 and n 3 / 2 , the spin decouples from the axis. Consequently, the ^-doubling 



is transformed into A-doubling. Equations (16) show that this can cause significant changes 



in sensitivity coefficients. The spin-orbit constant A can be either positive (CH molecule), 
or negative (OH). The sign of the fi-doubling depends on the sign of A, while A-doubling 
does not depend on A at all. Therefore, decoupling of the spin can change the sign of the 



splitting. In Sec. Ill B we will see that this can lead to a dramatic enhancement of the 



sensitivity to the variation of fundamental constants. 



B. Intermediate coupling 



The A-doubling for the intermediate coupling was studied in detail in many papers, 
including [63-65] (see also the book [59]). Here we use the effective Hamiltonian H s r from 



in the subspace of the levels n^ 2 and II3 , 2 , where upper sign corresponds to the parity p in 



Eq. (12). The operator H e g includes spin-rotational and hyperfine parts 

H c s = H BI + Hm ■ (17) 

Neglecting third order terms in the Coriolis and spin-orbit interactions, we get the following 
simplified form of the spin-rotational part: 

(n 1/2 , J,p|# sr |n 1/2 , J,p) = -\A + B(J + \f + p(S 1 + 5 2 )(2J + 1) , (18a) 

(U 3/2 ,J,p\H sr \U 3/2 ,J,p) = +lA + B(J+l) 2 -2B, (18b) 



(Tl 3/2 ,J,p\H SI \Tl 1/2 ,J,p)= [B + pS 2 {J+\)]^{J-\){J+l). (18c) 

Here in addition to the parameters A and B we have two parameters which appear in 
the second order of perturbation theory via intermediate state(s) Ei/ 2 . The parameter S\ 
corresponds to the cross term of the perturbation theory in the spin-orbit and Coriolis 
interactions, while the parameter S 2 is quadratic in the Coriolis interaction. Because of this 
S\ scales as o?[i and S 2 scales as /i 2 . It is easy to see that the Hamiltonian H sr describes 



limiting cases \A\ 3> B and \A\ <^ B considered in Sec. Ill A 



TABLE I: Frequencies (in MHz) and sensitivity coefficients for hyperfine components (J, F — >■ J, F') 
of A-doublet lines in CH and OH molecules. Recommended frequencies and their uncertainties are 
taken from I60H621 . 



Molecule Level J F F' to (MHz) Q a Q, 













Recom. 


Uncert. 


Theory 


Diff. 






12 CH 


2 n 1/2 


0.5 





1 


3263.795 


0.003 


3269.40 


-5.61 


0.59 


1.71 






0.5 


1 


1 


3335.481 


0.001 


3340.77 


-5.29 


0.62 


1.70 






0.5 


1 





3349.194 


0.003 


3354.11 


-4.92 


0.63 


1.69 






1.5 


1 


2 


7275.004 


0.001 


7262.25 


12.75 


-0.24 


2.12 






1.5 


1 


1 


7325.203 


0.001 


7312.02 


13.18 


-0.23 


2.11 






1.5 


2 


2 


7348.419 


0.001 


7335.30 


13.12 


-0.22 


2.11 






1.5 


2 


1 


7398.618 


0.001 


7385.08 


13.54 


-0.20 


2.10 


12 CH 


2 n 3/2 


1.5 


2 


2 


701.68 


0.01 


682.96 


18.72 


-8.44 


6.15 






1.5 


1 


2 


703.97 


0.03 


679.83 


24.14 


-8.66 


6.32 






1.5 


2 


1 


722.30 


0.03 


702.98 


19.52 


-8.37 


6.17 






1.5 


1 


1 


724.79 


0.01 


699.85 


24.94 


-8.07 


5.97 


16 OH 


n 3/2 


1.5 


1 


2 


1612.2310 


0.0002 


1595.42 


16.81 


-1.27 


2.61 






1.5 


1 


1 


1665.4018 


0.0002 


1648.93 


16.47 


-1.14 


2.55 






1.5 


2 


2 


1667.3590 


0.0002 


1650.66 


16.70 


-1.14 


2.55 






1.5 


2 


1 


1720.5300 


0.0002 


1704.17 


16.36 


-1.02 


2.49 


16 OH 


2 n 1/2 


0.5 





1 


4660.2420 


0.0030 


4638.98 


21.26 


2.98 


0.50 






0.5 


1 


1 


4750.6560 


0.0030 


4729.51 


21.15 


2.96 


0.51 






0.5 


1 





4765.5620 


0.0030 


4744.50 


21.06 


2.96 


0.51 






4.5 


5 


4 


88.9504 


0.0011 


64.34 


24.61 


-921.58 


459.86 






4.5 


5 


5 


117.1495 


0.0011 


92.35 


24.80 ■ 


-699.65 


349.59 






4.5 


4 


4 


164.7960 


0.0011 


141.20 


23.60 ■ 


-496.67 


248.77 






4.5 


4 


5 


192.9957 


0.0011 


169.22 


23.78 


-424.05 


212.68 



The hyperfine part of the effective Hamiltonian is defined in the lowest order of pertur- 
bation theory and has the form: 

(n 1/2 , J,p\H hf \U 1/2 , J,p) = C F [2a-b- c + p{2J + l)d] , (19a) 

(n 3/2 , J,p\H hi \U 3/2 , J,p) = 3C F [2a + b + c] , (19b) 

(n 3/2 , J,p|# hf |II 1/2 , J,p) = -C FV / (2J-l)(2J + 3)b, (19c) 

C F = [F(F + 1) - J{J + 1) - /(/ + 1)][8J( J + I)]' 1 . 

Here we assume that only one nucleus has spin and include only magnetic dipole hyperfine 
interaction. 

The effective Hamiltonian described by Eqs. ( 18|19 ) has 8 parameters. We use NIST 



values [SD] for the fine structure splitting A, rotational constant B, and magnetic hyperfine 
constants a, b, c, d. Remaining two parameters S\ and £2 are found by minimizing the rms 
deviation between theoretical and experimental A-doubling spectra. 
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FIG. 1: Sensitivity coefficients Q a and Q^ for A-doublet lines with AF = in CH and OH. The 
difference between lines with F = J + \ and F = J — ^ is too small to be seen. For the state n 3 / 2 
of OH the values for J = | are too large to be shown on the plot. They are listed in Table il 

In order to find sensitivity coefficients Q a we calculate transition frequencies for two 
values of a = a ± 5 near its physical value cto = 1/137.035999679(94). The similar 
procedure is applied to Q M at the physical value of the electron-to-proton mass ratio, 
/zo = 1/1836.15267247(80). We use scaling rules discussed above to recalculate parame- 
ters of the effective Hamiltonian for different values of fundamental constants. Then we use 
numerical differentiation to find respective sensitivity coefficient. 



C. Sensitivity coefficients for A-doublet transitions in CH and OH 



In Ref. [55J, the method described in the previous section was applied to 16 OH, 12 CH, 
7 Li 16 0, 14 N 16 0, and 15 N 16 0. The molecules CH and NO have ground state 2 II 1/2 (A > 0), 
while OH and LiO have ground state 2 n 3 / 2 (A < 0). The ratio |v4/-B| changes from 2 for CH 
molecule [UB], to 7 for OH [67J, and to almost a hundred for LiO and NO. Therefore, LiO 
and NO definitely belong to the coupling case a. For OH molecule we can expect transition 
from case a for lower rotational states to case b for higher ones. Finally, for CH we expect 
intermediate coupling for lower rotational states and coupling case b for higher states. 

Let us see how this scheme works in practice for the effective Hamiltonian ( l~8fl9 ). Fig. 
[l] demonstrates J-dependence of the sensitivity coefficients for CH and OH molecules. Both 
of them have only one nuclear spin I = |. For a given quantum number J, each A- 
doublet transition has four hyperfine components: two strong transitions with AF = and 
F = J ± | (for J = \ there is only one transition with F = 1) and two weaker transitions 
with AF = ±1. The hyperfine structure for OH and CH molecules is rather small and 
sensitivity coefficients for all hyperfine components are very close. Because of that Fig. [T] 
presents only averaged values for strong transitions with AF = 0. 

We see that for large values of J the sensitivity coefficients for both molecules approach 
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FIG. 2: Rotational spectrum of CH from (SB]. Vertical and diagonal arrows correspond to pure 
rotational and mixed transitions, respectively. A-doubling is not to scale. 



limit (16c) of the coupling case b. The opposite limits ( 16a|16b ) are not reached for either 
molecule even for smallest values of J. So, we conclude that the coupling case a is not 
realized. It is interesting that in Fig. [I] the curves for the lower states are smooth, while 
for upper states there are singularities. For CH molecule this singularity takes place for the 
state II3/2 near the lowest possible value J = 3/2. A singularity for OH molecule takes place 
for the state H±/2 near J = 9/2. 

These singularities appear because A-splitting turns to zero. As we saw above, the sign 
of the splitting for the coupling case a depends on the sign of the constant A. The same sign 
determines which state n^, or n 3 / 2 lies higher. As a result, for the lower state the sign of 
the splitting is the same for both limiting cases, but decoupling of the electron spin S for 
the upper state leads to the change of sign of the splitting. Of course, these singularities are 
most interesting for our purposes, as they lead to large sensitivity coefficients which strongly 
depend on the quantum numbers. Note, that when the frequency of the transition is small, 
it becomes sensitive to the hyperfine part of the Hamiltonian and the sensitivity coefficients 
for hyperfine components may differ significantly. The sensitivity coefficients of all hyperfine 
components of such A-lines are given in Table |T} We can see that near the singularities all 
sensitivity coefficients are enhanced. 

In addition to A-doublet transitions and purely rotational transitions there are also mixed 
transitions between rotational states of H1/2 and n 3 / 2 states. The transition energy here in- 
cludes the rotational and the fine structure parts. Because of that, such transitions may have 
different sensitivities to the variation of fundamental constants [56]. As an example, Fig.[2] 
shows mixed transitions in CH molecule. The sensitivity coefficients are given in Table |HJ 
The isotopologue CD has mixed transitions of lower frequencies and higher sensitivities [56J . 
Similar picture takes place for OH molecule. 

The molecule NH + is isoelectronic to CH and also has ground state 2 H.\/2- However, 
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TABLE II: Frequencies (GHz) and sensitivities of the rotational and mixed transitions in CH. 



N,J,p N',J',pf ftheor fcxpt [62J Q c 
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J-i 2 ' ' 
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1 ^ 
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1663.0 

1658.8 

2011.8 

2007.8 

193.1 

180.9 
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0.00 
0.00 
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0.42 
0.01 
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Ql 
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FIG. 3: Spin-rotational levels of the three lowest electronic states of the molecule NH + . Panels 
(a) and (b) correspond to vibrational states v = and v = 1 respectively. The energy levels are 
labeled with the quantum number J for the II states and with J and N for the £ state. 



there is an important difference: for NH + the first excited state 4 S~ lies only 340 cm -1 
above the ground state [681 EH]- The spin-orbit interaction between these states leads to 
strong perturbations of the rotational structure and of the A-doublet splittings and to an 
additional enhancement of the sensitivity coefficients [HE]. The spectrum of NH + is shown in 
Fig. |3j The effective Hamiltonian is similar to the one considered above with two additional 
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TABLE III: Low lying energy levels for the bending mode uj v = 589 cm -1 of Z-C3H molecule and 
their sensitivities q a and q^ to the variation of a and \i respectively. A is the distance from the 
ground state. All values are in cm -1 . 



"nom 


(v) K 


n 


(A) E 


A 
[73] 


% 


q a 





1.22 1 


0.5 


0.50 367.9 0.0 


0.0 


187.8 • 


-14.6 





1.35 1 


1.5 


0.46 381.9 13.9 


14.0 


187.8 


13.3 


1 


2.32 


0.5 


-0.01 394.2 26.3 


27.0 


197.3 


-0.4 


1 


3.57 2 


1.5 


0.21 597.7 229.7 


226.0 


300.3 


-6.1 


1 


3.65 2 


2.5 


0.19 603.5 235.5 


232.0 


300.3 


5.5 



terms describing interaction between the 2 I1 and 4 E states [68] : 

( 2 ir 3 /2, J,p|# so | 4 E~ 2 , J,p) = --C3/2 , (20a) 

( 2 n 1/2 , J,p|i7 so | 4 E- /2 , J,p) = -^Ci/2 • (20b) 

Obviously, the parameters C1/2 and C3/2 scale as a 2 . As mentioned above, for the NH + 
molecule the splitting between £ and II states AE^n is only about 340 cm -1 . This splitting 
includes three contributions: the non-relativistic electronic energy difference, the relativistic 
corrections (~ a 2 Z 2 ) and the difference in the zero point vibrational energies for the two 
states (~ /i 1 / 2 )- Note that the accidental degeneracy of these levels for NH + means that the 
first contribution is anomalously small. Because of that, the other two contributions can 
not be neglected and modify the scaling of AE^u with fundamental constants. This effect 
has to be taken into account in the calculations of the sensitivity coefficients [58] . 



IV. LINEAR POLYATOMIC RADICALS IN THE II GROUND STATE: C 3 H 

The linear form of the molecule C3H (Z-C3H) is similar to the molecule NH + : it also 
has the ground state 2 H±/2 and two closely lying states 2 n 3 / 2 and 2 St /2 . Here the quasi 
degeneracy of the II and £ states is not accidental, but is caused by the Renner- Teller 
interaction. In the following section we briefly recall the theory of the Renner- Teller effect 
in polyatomic linear molecules [70l [71] . 

A. Renner- Teller effect 

The total molecular angular momentum of the polyatomic molecule J includes the vi- 
brational angular momentum G associated with the twofold degenerate bending vibration 
mode(s): J = N+S = R+G+L + S, where R describes rotation of the molecule as a whole 
and is perpendicular to the molecular axis (. Other momenta have nonzero (^-projections: 
(G c ) = I, (L c ) =A,(N c ) = K = l + A, and ( J c ) = Q. 

Suppose we have II electronic state |A = ±1) and v = 1 vibrational state of a bending 
mode |Z = ±1). All together there are 4 states |A = ±1)|Z = ±1). We can rewrite them 
as one doublet A state \K = ±2) and states E + and E - . In the adiabatic approximation 
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all four states are degenerate. Renner [70] showed that the states with the same quantum 
number K = I + A strongly interact, so the S + and S~ states repel each other, while the A 
doublet in the first approximation remains unperturbed. We are particularly interested in 
the case when one of the £ levels is pushed close to the ground state v = 0. This is what 
takes place in the /-C 3 H molecule [TIHTS] . 

Consider a linear polyatomic molecule with the unpaired electron in the 7r^ state in the 
molecular frame £, r), (. Obviously, the bending energy is different for bendings in £( and in 
rj( planes: V± = \k±\ 2 (here \ is the supplement to the bond angle). That means that the 
electronic energy depends on the angle (f> between the electron and nuclear planes: 



H' = V'cos2(j), 

where 2V = V + — V- — k'x 2 ■ There is no reason for V to be small, so k' ~ k^ 
to a first approximation k! does not depend on a and \i. 



(21) 
1 a.u. and 



As long as interaction (21) depends on the relative angle between the electron and the 
vibrational planes, it changes the angular quantum numbers as follows: AA = —A/ = ±2 
and AK = 0. This is exactly what is required to produce splitting between the S + and £~ 
states with v = 1 as discussed above. 



TABLE IV: Z-C3H sensitivity coefficients for the transitions between states from Table III and for 
parameters A e Q and A-E^n defined by (22) and (29) respectively. Frequencies are in cm -1 . 



k n K' Q' 



Fit to [73] 

w Qii Qc 



uj 



Fit to [76] 



1 


0.5 1 


1.5 13.9 


0.00 


2.00 14.4 


0.00 


2.00 


1 


1.5 


0.5 12.4 


0.78 


-1.11 13.3 


0.77 


-1.07 





0.5 2 


1.5 203.5 


0.51 


-0.03 204.4 


0.51 


-0.03 


2 


1.5 2 


2.5 5.8 


0.00 


2.00 6.0 


0.00 


2.00 




As 


13.9 


0.00 


2.00 14.4 


0.00 


2.00 




A-Esn 


19.4 


0.50 


0.00 20.5 


0.50 


0.00 



Interaction (21) also mixes different vibrational levels with Av = ±2, ±4, 



. Thus, 

we have, for example, the nonzero ME (0, 0, 1, 1\H' |2, 2, — 1, 1) between states \v,l,A,K). 
Such mixings reduce effective value of the quantum number A and, therefore, reduce the 
spin-orbital splitting between the II1/2 and n 3 / 2 states [73] , 



H so = A cff AS , A cS = AA cS /A . 



(22) 



Let us define the model more accurately. Following 

H = H e + T v + AL C S C 



we write the Hamiltonian as: 

(23) 



Here the "electronic" part H e includes all degrees of freedom except for the bending vibra- 
tional mode and spin. For Z-C3H there are two bending modes, but for simplicity we include 
the second bending mode in H e too. Electronic MEs in the |A) basis have the form: 



(±l|# e |±l) 
(±l|ff e |=Fl) 



V+ + V- k 



;X 



k' 



X 2 exp (=p2i<; 



(24a) 
(24b) 
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Here x an d 4> are the vibrational coordinates for the bending mode. Kinetic energy in these 
coordinates has the form: 

2MR 2 \d X 2 xdxX 2 d(j) 2 )' { ' 

We can use the basis set of 2D harmonic functions in polar coordinates p = xR an d <fr 
for the mass M and the force constant k: 

Vw (P, 0) = Rv,i (p) —?= exp (il(j)) . (26) 

v 2ir 



It is important that the radial functions are orthogonal only for the same I: 

(Rv',l\Rv,l) — $v',v ■ (27) 

This allows for the nonzero MEs between states with different quantum number /. By 



averaging operator (23) over vibrational functions we get: 



(v',l'\H e + T v \v,l) = [oj v (v + 1)+ AAS^]S V ' ;V S, 



+ ^(R v >i'\k'x 2 \Rvi) exp (t2Wm±2 • (28) 

The exponent here ensures the selection rule A' = A =p 2 for the quantum number A when 
we calculate MEs for the rotating molecule. 

B. Molecule Z-C,H 



We solve the eigenvalue problem for Hamiltonian (23) using the basis set of the 2D- 
harmonic oscillator. Our model Hamiltonian has only 3 parameters, namely u v , A, and 
the dimensionless Renner- Teller parameter S: k! = Sk. The values for cj v and A for I-C3H. 
are given in [73]. We varied the Renner- Teller parameter £ to fit five lowest levels for the 
given bending mode: n^, n 3 / 2 , £1/2, A 3 / 2 , and A 5 / 2 . The optimal value appeared to be 



£ = 0.788. The results are presented in Table III The first two columns give nominal 



vibrational quantum number v and its actual average value. We see that the Renner- Teller 



term in (28) strongly mixes vibrational states. This mixing also affects (A) and decreases 



spin-orbital splittings as explained by Eq. (22). 



cm 



The last two columns in Table III give dimensional sensitivity coefficients q^ and q a in 

1. — 



A rn AC* AU 

AE = q a + <^— . 

a jjl 

To estimate them we assumed that the parameters scale in a following way: oj v ~ \x 



1/2 



A ~ a , and £ does not depend on a and /1. The dimensionless sensitivity coefficients (10) 
for the transitions w^ = E\. — Ei can be found as: 

Qi,k = {qk — Qi)/^i,k ■ 



In Table [TV] these coefficients are calculated for the same set of parameters as in Table [ITT 
and for the slightly different parameters which better fit experimental frequencies from [76J . 
We see that the sensitivity coefficients are practically the same for both sets. 
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FIG. 4: Spin-rotational levels of the three lowest vibronic states of the molecule Z-C 3 H. i^-doubling 
is indicated schematically, hyperfine structure is shown only for the two lowest -fC-doublets. Due 
to a strong Renner- Teller effect the component 2 £ + of the excited bending state 1/4 (CCH bending) 
is shifted towards lower energies, ~ 29 cm -1 above the zero-level of the ground state 2 n]/ 2 . 



For the two fine structure transitions, IIi/ 2 — > n 3 / 2 and A 3 / 2 — > A 5 / 2 , we get sensitivi- 
ties Qn = and Q a = 2. This may seem strange as the fine structure is significantly reduced 
by the Renner- Teller mixing: the fine-structure parameter is 29 cm -1 and the splitting be- 
tween II1/2 and n 3 / 2 is only 13.9 cm" 



According to (22) the mixing reduces the splitting. 



However, this effect depends on the dimensionless Renner- Teller parameter £ and does not 
depend on \i and a. Consequently, the effective parameter A e s depends on fundamental 
constants in the same way as initial parameter A. 



For the high frequency transition E]/ 2 



A3/2, where the spin-orbital energy can be 



neglected, we get Q^ = 0.5 and Q a = 0. These results are expected, because our model has 
only two dimensional parameters: vibrational frequency, which is proportional to /i 1//2 and 
the fine structure parameter A, which scales as a 2 . Even though our vibrational spectrum is 



far from that of a simple harmonic oscillator, the non-diagonal MEs (28) of the Hamiltonian 
(23) still scale as /i 1 / 2 . Therefore, if we neglect spin-orbital splittings, we get Q^ = 1/2 for all 



transitions. The only transition in Table IV where the spin-orbital energy and vibrational 
energy are close to each other is the n 3 / 2 — > Ex/ 2 transition. The resultant frequency is 
roughly half of the vibrational energy difference between the II and S states. This leads to 
Qf, rs 1 and Q a fa -1. 

The spectrum of the /-C 3 H molecule is shown on Fig. |4j The effective Hamiltonian 
for the rotating molecule is similar to that of the NH + molecule. It includes the effective 
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TABLE V: Frequencies (MHz), sensitivity coefficients, and reduced MEs (a.u.) for some K- 
doubling transitions in n 1 / 2 state of the Z-C3H molecule. 



J F'p, Fp 



UJ 



Q, 



Q f 



\D\ 



il+, 


0- 


52.37 


0.66(2) 


1.7(2) 


0.333 


|o+, 


1- 


39.12 


0.20(2) 


1.9(2) 


0.333 


£l+, 


1- 


34.93 


-0.02(2) 


2.0(2) 


0.667 


§1-. 


1+ 


85.55 


0.65(2) 


1.7(1) 


0.166 


5 2- 

2 z ' 


1+ 


78.60 


0.55(2) 


1.7(1) 


0.033 


2 L > 


2+ 


75.23 


0.43(2) 


1.8(1) 


0.033 


5 2- 

2 z ' 


2+ 


68.29 


0.30(2) 


1.8(1) 


0.299 


|2+, 


2- 


107.19 


0.95(2) 


1.5(1) 


0.132 


|3+, 


2- 


98.97 


0.89(2) 


1.5(1) 


0.009 


|2+, 


3- 


98.83 


0.82(2) 


1.6(1) 


0.009 


|3+, 


3- 


90.61 


0.75(2) 


1.6(1) 


0.188 


2 J ' 


3+ 


112.38 


1.63(2) 


1.2(1) 


0.105 


14-, 


4+ 


96.07 


1.56(2) 


1.2(1) 


0.136 


§4+, 


4- 


95.75 


3.22(4) 


0.36(7) 


0.086 


15+, 


5- 


79.63 


3.45(4) 


0.23(7) 


0.105 


11 5- 
2 D 


, 5+ 


52.81 


9.1(6) 


-2.6(3) 


0.072 


^6- 
2 u 


, 6+ 


36.85 


12.1(6) 


-4.1(3) 


0.085 


2 u 


, 6+ 


20.25 


-34.(2) 


19.(2) 


0.062 


13 7 
2 ' 


, 7+ 


36.06 


-18.(2) 


11.(2) 


0.071 


f7+ 


, 7- 


126.59 


-7.6(2) 


5.8(4) 


0.054 


f8+ 


, 8- 


142.24 


-6.5(2) 


5.3(4) 


0.061 


¥«- 


, 8+ 


268.76 


-4.7(1) 


4.4(3) 


0.047 


¥9" 


, 9+ 


284.25 


-4.3(1) 


4.2(3) 


0.053 


f9+ 


, 9- 


448.75 


-3.59(7) 


3.8(3) 


0.042 


fio+ 


, 10- 


464.07 


-3.39(7) 


3.7(3) 


0.046 


fio- 


, 10+ 668.02 


-2.97(6) 


3.5(3) 


0.038 


fll- 


, 11+ 683.18 


-2.85(6) 


3.4(3) 


0.041 



fine-structure parameter A e g and the energy difference between the £ and II states, 



AE- 



ai 



£(E" 



£(II 1/2 ) + E{U 3/2 ) 



(29) 



Numerical values for these parameter are obtained from the fit to the experimental transi- 
tion frequencies. Here we onl y n eed to determine the dependence of these parameters on 



fundamental constants. Table IV shows that A c g ~ a and AE^u ~ jjl ' . Once again, this 



is because the Renner- Teller mixing depends on the dimensionless parameter £ and does not 
depend on a and fi. Calculated sensitivity coefficients for the i^-doublet transitions of the 



Z-C3H molecule are listed in Tables [V] and |VI[ The results for the mixed transitions can be 
found in 
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TABLE VI: Frequencies (MHz), sensitivity coefficients, and reduced MEs (a.u.) for some K- 
doubling transitions in II3 / 2 state of the Z-C3H molecule. 



JF'p',Fp 



u> 



^ca 



Q> 



\D\ 



§1-1+ 


5.61 


-2.63(8 


3.2(2) 


1.493 
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0.299 
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22.24 
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§3+,2- 


31.50 


-1.35(8; 
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12.88 
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0.085 
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22.15 
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1.694 
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54.92 
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14-, 4+ 


54.76 
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5 +- 
2 ^ 


108.13 


-2.50(8 


3.1(2) 


1.230 


11 -+ 
2 t 


185.99 


-2.46(8; 


3.1(2) 


1.007 


39 1 

2 " l " 


4266.17 


-2.9(1) 


2.53(8) 


0.224 


^+- 
2 ~ 


4553.04 


-3.5(1) 


2.42(5) 


0.208 


43 _, 
2 ' 


4663.43 


-4.6(2) 


2.2(1) 


0.192 


45 ,_ 
2 ~ 


4377.16 


-7.5(2) 


1.4(3) 


0.174 


47 _, 
2 T 


3097.96 


-19.0(4) 


-2.3(9) 


0.149 


49 1 
2 "i" 


909.06 


132.(2) 


53.(8) 


0.103 


51 _ + 

2 ~ 


19813.69 


-3.11(5 


-1.6(4) 


0.116 



V. TUNNELING MODES IN POLYATOMIC MOLECULES 



In this section we consider non linear and non planar polyatomic molecules. Such 
molecules generally have more than one equivalent potential minimum. If the barriers be- 
tween these minima are not too high the molecule can tunnel between them. Ammonia (NH3) 
is the best known textbook example of a nonrigid molecule (see Fig. pjj). Interestingly, this 
molecule is also one of the most abundoned polyatomic molecules in the interstellar medium. 
Other important for astrophysics molecules with tunneling include hydronium (H 3 + ), per- 
oxide (H2O2), methanol (CH3OH), and methylamine (CH3NH2). We will briefly discuss all 
of them below. All these molecules include only light atoms with Z < 8 and have singlet 
electronic ground states. Thus we can neglect relativistic corrections and assume that all 
discussed transitions have Q a = 0. 

It is clear that tunneling frequencies should strongly depend on the nuclear masses, and 
we can expect large sensitivity coefficients Q^tw.- They can be found using the semi-classical 
Wentzel-Kramers-Brillouin (WKB) approximation. Following [77] we can write the ground 
state tunneling frequency in atomic units (h — \e\ = m e = 1) as: 



2E, 



Wu- 



-s 



7T 



(30) 



where S is the action over classically forbidden region and E is the ground state vibrational 
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FIG. 5: Potential for the tunneling (umbrella) mode of the NH3 molecule. Two lowest vibrational 
levels lie below the barrier. 

energy calculated from the bottom of the well U m \ n . If the barrier is high enough the 
harmonic approximation gives 2Eq = u v , where u v is the observed vibrational frequency. In 
this case Eq. (30) allows to find action S from experimentally known frequencies a; tun and u v . 
For lower barriers we need to know the shape of the potential to estimate Eq. The examples 
of these two limiting cases are ammonia and hydronium, where tunneling frequencies are 0.8 
cm -1 and 55 cm -1 respectively. 



/' 



The action 5* depends on the tunneling mass, which in atomic units is proportional to 
_1 . Differentiating (30) over /i we get [751 179"] 



Q 



s 



SE 



/i,tun 



2 (AC/ - E 



(31) 



where AC/ = C/ max — C/ m i n is the barrier hight. Numerical solution of the Schrodinger equation 
for realistic potentials agrees with this WKB expression within few percent for all molecules 
considered so far. 



A. Ammonia 



Equations (30), (31) show that sensitivity coefficient logarithmically depends on the tun- 



neling frequency. For example, for the symmetric isotopologues of ammonia we get: 



NH 3 : 
ND 3 : 



w tun = 24 GHz, 
u tun = 1.6 GHz, 



W/x,tun = u. l • 



(32a) 
(32b) 



Such a weak dependence on the tunneling frequency limits possible values of the sensitivity 
coefficients for tunneling transitions in the microwave range: Q^tun ^ 8. This is quite good, 
compared to the rotational sensitivity Q^rot ~ 1, but smaller than the best sensitivities in 
linear molecules considered above. 

Let us consider mixed tunneling-rotational transitions, where tunneling goes along with 
the change of the rotational quantum numbers. If we neglect interaction between tunneling 



19 



and rotational degrees of freedom we can write approximate expressions for the frequency 
and the sensitivity of the mixed inversion-rotational transition: 

^rnix = Wri ^tun , (33a) 

k'mix ^mix 



We are particularly interested in the case when the minus sign in (33) is realized and 
^mix *C k>tun- For this case the tunneling sensitivity is enhanced by the factor (jJ tun /uj mix ^> 1 
and resultant sensitivity of the mixed transition is inversely proportional to the transition 
frequency u m i x . Therefore, for the mixed transitions we can have much higher sensitivities 
in the observable frequency range, then for the purely tunneling transitions. 

Another important advantage of the mixed transitions is that there are usually many of 
them each having different sensitivity. This means that we can have very good control on 
possible systematics and reliably estimate the accuracy of the results for /i-variation. 

The mixed transitions can not be observed in the symmetric isotopologues of ammonia 



(32 ), but they are observed in the partly deuterated species NH 2 D and NHD 2 . Unfortunately, 



for both of them the tunneling frequency is much smaller then all rotational frequencies and 



sensitivities (|33b| are not large 

(34) 



NH 2 D: 0.10 <g M , mix < 1.61, 



NHD 2 : 0.27 < Q„ mix < 1.54 . 



Relatively small sensitivity coefficients for deuterated isotopologues of ammonia (34) and 



their low abundance does not allow to get strong limits on /i-variation, so we need to use 



tunneling ammonia line (32a). It was observed from the several objects with the redshifts 
about unity. Measuring radial velocities for rotational lines and for the ammonia tunneling 
line we have AQ = 3.5 in Eq. (pi), which is two orders of magnitude larger than for optical 
lines. Because of that the ammonia method allowed to place more stringent bounds on 
/^-variation than bounds, which follow from the optical spectra of the hydrogen molecule. 
However, recent observations of the molecules with mixed tunneling-rotational transitions 
provide even higher sensitivity to /i-variation. 

B. Mixed tunneling-rotational transitions and effective Hamiltonians 



Equations (33) show that high sensitivity mixed transitions are possible when tunneling 
frequency is of the same order of magnitude as rotational constants. However, in this case 
tunneling and rotational degrees of freedom start to interact and the accuracy of approxi- 



mation (33) decreases. A much better approximation can be reached with the help of the 
effective Hamiltonians, which describe rotational and tunneling degrees of freedom and their 
interactions with each other. At present the state of the art effective Hamiltonians can in- 
clude on the order of hundred parameters. These parameters are fitted to the experimentally 
known transitions and provide an accuracy on the ppm scale, or better. 

When such Hamiltonians are used to find sensitivity coefficients Q^ we need to know 
how all the parameters depend on /i. It was shown in [25] that this can usually be done 
only within an accuracy of a few percent. The final accuracy for the large Q" factors i s 



somewhat lower because of the instability of Eq. (33b ). Because of that we need not complex 
effective Hamiltonians but their simplified versions with considerably smaller numbers of 
fitting parameters can be used instead. 
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FIG. 6: Tunneling-rotational spectrum of HsO + molecule. Several low-frequency tunneling and 
mixed transitions are marked with vertical arrows. Their frequencies are shown in GHz. 

C. Hydronium and peroxide 



Let us start with hydronium molecule H3O 4 " [HI]. This molecule is a symmetric top. It is 
similar to ammonia, but natter. Tunneling frequency is almost 50 times larger and compa- 
rable to rotational intervals. The tunneling umbrella mode does not change the symmetry 
and does not contribute to the angular momentum of the molecule. Because of that the 
tunneling-rotational interaction is reduced to the centrifugal corrections to the tunneling 
frequency [82] . 

The tunneling-rotational spectrum of hydronium is shown in Fig. [6] It consists of the 
J ladders for each quantum number K, where K is projection of the angular momentum 
on the molecular axis. Due to the tunneling each rotational level is split in two states with 
different parity p. For K = the permutation symmetry of the hydrogen nuclei allows only 
one of these levels, while for K > both levels are present. 

In Fig. [6] we see four mixed transitions with frequencies around 300 GHz, which is few 
times smaller than the tunneling frequency that is about 1.6 THz. Table |VTT shows that these 
transitions have enhanced sensitivity to \i variation (Q^tun — 2.0 ±0.1). Those transitions, 
whose frequencies decrease when tunneling frequency increases have negative sensitivity 
coefficients Q^. We conclude that hydronium has several mixed transitions with sensitivities 
of both signs and the maximum AQ^ is around 10. Other isotopologues of hydronium have 
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TABLE VII: Sensitivities of the low frequency mixed inversion-rotational transitions in hydronium 
H30 + . Molecular states are labeled with quantum numbers J V K . 



Transition 


Frequency (MHz) 


Q» 


Upper Lower 


Theory Exper. 




If 


n 


307072 307192.4 


6.4(5) 


3+ 


^ 


365046 364797.4 


-3.5(5) 


n 


2r 


389160 388458.6 


-3.1(4) 


°0 


2n 


397198 396272.4 


-3.0(4) 


On 


b! 


984690 984711.9 


2.7(2) 


4 


3* 


1031664 1031293.7 


-0.6(2) 


4 


3-7 


1071154 1069826.6 


-0.5(2) 


32 


3^ 


1621326 1621739.0 


2.0(1) 


2f 


2| 


1631880 1632091.0 


2.0(1) 


ir 


If 


1655832 1655833.9 


2.0(1) 



even higher sensitivities [79], but up to now they have not been observed in the interstellar 
medium. 

Another molecule where tunneling frequency is comparable to rotational constants, but 
tunneling-rotational interaction is rather weak, is peroxide H2O2 [E31 EI]- In equilibrium 
geometry H 2 2 is not flat; the angle 27 between two HOO planes is close to 113°. Two flat 
configurations correspond to local maxima of potential energy; the potential barrier for trans 
configuration (27 = it) is significantly lower, than for cis configuration (7 = 0): U n ~ 400 
cm -1 and Uq ~ 2500 cm -1 . To a first approximation one can neglect the tunneling through 
the higher barrier. In this model peroxide is described by a slightly asymmetric oblate top 
with inversion tunneling mode, similar to ammonia and hydronium. 

The sensitivity coefficients for the mixed transitions in peroxide were calculated in [86] . 
Results of these calculations are shown in Table IVIIII Molecular states are labeled with 
the rotational quantum numbers J, K&, and Kq and the tunneling quantum number r [83J. 
Transitions with the frequencies below 100 GHz were found to have rather high sensitivities of 
both signs. Several transitions of peroxide were recently observed from interstellar medium 
(ISM) in [85]. These transitions have higher frequencies and smaller sensitivities to fi- 
variation. Nevertheless, even for these transitions the maximum value of AQ^ is about 
4.5. 



D. Molecules with hindered rotation: methanol and methylamine 

Hindered rotation is one of the examples of the large amplitude internal motions in non 
rigid molecules. In the discussion of the peroxide molecule in the previous subsection, we 
neglected the tunneling through the higher cis barrier. For the excited vibrational states 
tunneling through both barriers can take place leading to the hindered rotation of one 
HO group in respect to another. Many molecules which include CH 3 group have three 
equivalent minima at 120° to each other. Hindered rotation in such molecules can take place 
already for the ground vibrational state. When the tunneling frequencies are comparable to 
the rotational ones, such molecules have very rich microwave spectra with a large number 
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TABLE VIII: Numerical calculation of the Q-factors for low frequency mixed transitions in peroxide 
H2O2 using effective Hamiltonian. Experimental frequencies are taken from JPL Catalogue |61| . 
E up is upper state energy in Kelvin. 



Jk a ,k c (t) 


-E'up 


UJ 


(MHz) 


Q» 


upper 


lower 


(K) 


theory 


exper. 






Transitions below 100 GHz 




0o,o(3) - 


-li,o(l) 


17 


14818.8 


14829.1 


+36.5(2.9) 


214(1) - 


- lo,i(3) 


21 


37537.0 


37518.28 


-13.0(1.2) 


lo,i(3) - 


-li,i(l) 


19 


67234.5 


67245.7 


+8.8(6) 


2o, 2 (3) - 


-2i )2 (l) 


24 


68365.3 


68385.0 


+8.7(6) 


3o, 3 (3) - 


-3i )3 (l) 


31 


70057.4 


70090.2 


+8.5(6) 


4o,4(3) - 


- 4l,4(l) 


41 


72306.0 


72356.4 


+8.3(6) 


5o, 5 (3) - 


- 5i, 5 (l) 


53 


75104.6 


75177.4 


+8.0(6) 


6o, 6 (3) - 


-6i,6(l) 


68 


78444.7 


78545.4 


+7.7(6) 


3i, 2 (l) - 


- 2 0i2 (3) 


28 


90399.8 


90365.51 


-4.8(5) 




Transitions observed from ISM in 


[85] 


3o,3(3) - 


-2i,i(l) 


31 


219163.2 219166.9 


+3.4(2) 


6i,5(l) - 


- 5o, 5 (3) 


66 


252063.6 251914.68 


-1.1(2) 


4o,4(3) - 


-3i, 2 (l) 


41 


268963.7 268961.2 


+3.0(2) 


5o, 5 (3) - 


- 4i )3 (l) 


53 


318237.7 318222.5 


+2.7(1) 



of mixed transitions. Another distinctive feature of these molecules is strong interaction 
between the internal (hindered) and overall rotations. One of the simplest molecules of this 
type is methanol CH3OH. 

The basic theory of the non-rigid tops with internal rotation was established in the 
1950s JH3 EH] and the main features of the methanol spectrum were explained. Later on the 
theory was refined many times and currently there is a very impressive agreement between 
the theory and experiment [89H92] . 

The sensitivity coefficients to the /^-variation for methanol microwave transitions were 
calculated independently in [531 El] and in [2S]. The first group used the state of the art 
effective Hamiltonian |91j . which included 120 fitting parameters. The second group used a 
much simpler model [HS] ■ The rotational part H mt was that of the slightly asymmetric top 
and included the rotational constants A, B, and C (A « B). The hindered rotation was 
described by the Hamiltonian 



H\ 



hr 



r, d2 V 3 

-F h — 

duo 2 2 



cos 3w) 



(35) 



where the kinetic coefficient F was proportional to \i and the electronic potential V3 was 
independent on \i. The angle uj described position of the OH group in respect to the CH 3 top. 
This model did not include centrifugal distortions. The interaction of the internal rotation 
with the overall rotation was described by a single parameter D, which scaled linearly with 
\i [HZ]. Altogether this model had 6 parameters. 

Both effective Hamiltonians were diagonalized for several sets of parameters, which cor- 
respond to an increased and decreased // and the sensitivity coefficients were found by the 
numerical differentiation. The comparison of the two calculations is given in Fig. [7| We see 
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FIG. 7: Comparison of the sensitivity coefficients for CH3OH from [93] and [25]. The former used 
the sensitivity coefficients K^ defined as K^ = — Q^. This corresponds to the different definition 
of the mass ratio: m p /m e instead of m e /m p , which is used in the present review. 



that in spite of a significant difference in complexity of the models the results are in good 
agreement and of the comparable accuracy. As we discussed above, the latter is mostly 
determined by ambiguity in the /x-scaling of model parameters. 

The sensitivity coefficients for the mixed transitions in methanol span from —17 to +43, 
which corresponds to |AQ^| ~ 60. This is more than an order of magnitude larger than 
in ammonia method. Moreover, in methanol we have a large number of strong lines with 
different sensitivities and can effectively control possible systematic effects. Until very re- 
cently methanol was observed only at small redshifts, but in 2011 it was first detected in the 
microwave survey towards the object PKS 1830-211 at redshift z = 0.89 [96]. This means 
that at present methanol can be used as a very sensitive tool to probe /i-variation on a 
cosmological timescale [231 EZ] • 

In the same survey [96J, a large number of rather complex molecules were detected for 
the first time at high redshift. In particular, the list includes methylamine — yet another 
molecule with tunneling motion. In contrast to all previously discussed molecules, methy- 
lamine has two tunneling modes. First is hindered rotation of the NH2 around CH3 top, 
which is similar to that in methanol. Second is a wagging mode when the NH2 group flips 
over to the other side (see Fig. [8]). Both modes contribute to the angular momentum of the 
molecule and, therefore, strongly interact with the overall rotation. 

The spectrum of methylamine is also very rich. The effective Hamiltonian must include 
both tunneling motions and their interactions with each other and with the overall rotation. 
Therefore, even the simplest form of this Hamiltonian is quite complex and we will not 
discuss it here. Calculations of the sensitivity coefficients were recently done in [98]. It was 
found that they lie in the range —24 < Q M < 19. However, the lines, which were observed 
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FIG. 8: Schematic representation of methylamine and variation of the potential energy of methy- 
lamine as function of the relative rotation 7 of the CH3 group with respect to the amine group 
about the CN bond and the angle r of the two hydrogen atoms of the NH2 group with respect 
to the CN bond. The two large amplitude motions, corresponding to inversion hi v and hindered 
rotation h^ v are schematically indicated by the arrows. Note that inversion of the NH2 group is 
accompanied by a tt/3 rotation about the CN bond of the CH3 group with respect to the amine 
group. 



in [96] at z = 0.89 have sensitivities close to 1. Up to now neither of the more sensitive lines 
of methylamine has been observed at high redshifts. 

There are several other molecules with mixed tunneling-rotational spectra, for example 
N 2 H 4 and CH 3 SH. The former has three tunneling modes which strongly interact with 
rotation. Thus, we should expect very complex spectrum. This molecule is predicted to form 
in Jupiter's and Titan's atmospheres [9911100] . The latter is similar to methanol and exhibits 
hindered rotation and complex spectrum [101] . Effective Hamiltonians for many of these 



molecules are known, but no other calculations of the sensitivity coefficients have been done 
so far (preliminary results for CH3SH show that there are transitions with high sensitivities 
of both signs). If any new sufficiently low frequency mixed transitions are observed from 
the interstellar medium, it is possible to calculate respective sensitivity coefficients using the 
methods outlined in this section. 
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VI. SUMMARY AND CONCLUSIONS 

As we discussed in the previous sections the constraints on the possible variation of 
fundamental constants are an efficient method of testing the equivalence principle which is a 
basic assumption of General Relativity. These constraints can be derived from a wide variety 
of atomic and molecular transitions observed in laboratory, solar and extra solar systems, 
and at very early cosmological epochs up to a redshift of order z ~ 5 — 6 where molecular 
and atomic transitions have been recently detected and observed with a sufficiently high 
spectral resolution [TJ1 [15]. Radio astronomical observations of the NH3 molecule in two 
distant galaxies provide tight constraints at the A/x//z < 1 x 10~ 6 level at z = 0.89 [102] and 
z = 0.69 [13]. Even deeper bounds were deduced from observations of the CH 3 OH molecule 
in the z = 0.89 galaxy: A///// < 3 x 10~ 7 [97], and A///// < 1 x 10" 7 [23]. 

To probe a and \x at the level of 10 -8 or 10~ 9 , at least two main requirements should be 
fulfilled: {%) increasing precision of the laboratory measurements of the rest frame frequencies 
of the most sensitive molecular transitions discussed in this review, and {%%) increasing 
sensitivity and spectral resolution of astronomical observations. 

The most promising molecular transitions are those of a mixed nature, where there are 
two, or more, competing contributions to the transition energy. We get strong enhancement 
of the sensitivity to the variation of the fundamental constants when the resultant transition 
frequency is much smaller than individual contributions. This happens, for example, for 
some mixed tunneling-rotational transitions. Diatomic radicals give another example, where 
spin-orbit interaction is competing with Coriolis interaction. As a result we have strong 
enhancement of the sensitivity coefficients for the A-doublet transitions. There are other 
known examples, which are more relevant for the laboratory experiments [1031 1104] . It is 
possible that more examples will be found both for the laboratory and astrophysical studies. 
The methods described in this review allow us to calculate sensitivity coefficients for any 
microwave and submillimeter molecular transitions of interest. 
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